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Abstract 

We discuss recent developments in neutrino physics and focus, in particular, on 
neutrino oscillations and matter effects of three light active neutrinos. Moreover, 
we discuss the difference between Dirac and Majorana neutrinos, neutrinoless (5(5- 
decay, absolute neutrino masses and electromagnetic moments. Basic mechanisms 
and a few models for neutrino masses and mixing are also presented. 
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1 Motivation 



In recent years neutrino physics has gone through a spectacular development (for reviews 
see, for instance, PU 121 El 13])- Data concerning solar and atmospheric neutrino deficits have 
been accumulated, these deficits have been established as neutrino physics phenomena 
and the Solar Standard Model [5] has been confirmed. The last steps of this exciting 
development were the results of the SNO [S| and KamLAND experiments [7j: The SNO 
experiment provided a model-independent proof of solar v e — > v^ T transitions and the 
terrestrial disappearance of v e reactor neutrinos in the KamLAND experiment has shown 
that the puzzle of the solar neutrino deficit is solved by neutrino oscillations [Hj. This 
gives us confidence that the same is true for the atmospheric neutrino deficit as well. For 
information on the history of neutrino oscillations see j^l EI] , for the recent experimental 
history see the contribution of G. Drexlin to these proceedings [3]. General reviews on 
neutrino physics can be found, for instance, in Refs. jTT| H21 E3 EH] • 

These lecture notes are motivated by this development and aim at supplying the 
theoretical background for understanding and assessing it. In view of the importance 
of neutrino oscillations in this context, we will give a thorough discussion of vacuum 
neutrino oscillations and matter effects [T51 EH] (see Section |2J); the description of the 
latter is tailored for an understanding of the flavour transformation of solar neutrinos and 
effects in earth matter. Then, in Section El we will switch to the subject of the neutrino 
nature, which is a question of principal interest but has no impact on neutrino oscillations; 
we will work out the difference between Dirac and Majorana neutrinos and the basics of 
Majorana neutrino effects. Eventually, in Section 0] we will come to the least established 
field of neutrino physics: models for neutrino masses and mixing. In view of the huge 
number of models and textures, we cannot try to cover the field but rather discuss a small 
selection of basic mechanisms for generating neutrino masses and mixing. This selection 
will necessarily be biased due to personal interest and prejudices. A similar judgement 
has to be made concerning the literature quoted in these lecture notes; owing to the host 
of papers which have appeared in recent years only a small selection can be quoted here 
(for literature on neutrino experiments see Ref. Finally, we present conclusions in 
Section EJ 

Abbreviations used in these lecture notes: CC = charged current, NC = neutral cur- 
rent, LBL = long baseline, SBL = short baseline, MSW = Mikheyev-Smirnov-Wolfenstein, 
LMA = large mixing angle, MM = magnetic moment, EDM = electric dipole moment, 
SM = Standard Model, SUSY = supersymmetry, GUT = Grand Unified Theory, VEV = 
vacuum expectation value. 

2 Neutrino oscillations 

2.1 Neutrino oscillations in vacuum 

Here we give a simple and yet quite physical derivation of the neutrino oscillation formula. 
The first observation is that — as the NC interactions are flavour-blind — neutrino-flavour 
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production and detection proceeds solely via CC interactions with the Hamiltonian density 

n cc = ^=w; J2 LY^l + h.c. (i) 

The next observation is that flavour transitions are induced by neutrino mixing: 

VaL = J2 U »j U jL- (2) 

j 

This relation means that the left-handed flavour fields are not identical with the left- 
handed components of the neutrino mass eigenfields Uj corresponding to the mass rrij, 
but are related via a matrix U, which is determined by the neutrino mass term which will 
be discussed in Section 3. Here we simply assume the existence of the mixing matrix U 
and the neutrino masses rrij. 

We confine ourselves, apart from a few side remarks, to the following basic assump- 
tions: 

• There are three active neutrino flavours; 

• The mixing matrix U is a 3 x 3 unitary matrix. 

These assumptions are supported by the results of the neutrino oscillation experiments 

m 

A crucial observation is that neutrino flavour is defined by the associated charged 
lepton in production and detection processes; there is no other physical way to define 
neutrino flavours. Looking at Eqs. (JTJ and (J2J), we find that, in a reasonable approxima- 
tion, neutrino-flavour states are given by 



^> = E^iK>> ( 3 ) 
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y. 

3=1 

where the mass eigenstates \rrij) fulfill the orthogonality condition 

{m j \m k ) = 5 jk . (4) 

The state (j3J) describes a neutrino with flavour a sitting at the position x = 0. Every 
mass eigenstate propagates as a plane wave. If we have a stationary neutrino state with 
neutrino energy E, we derive from Eq. (jSJ) the following form of the propagating state: 

K,x) = j:U* aj e- i ( Et -^\m j }, (5) 



with neutrino momenta 



m 2 



Pj = y E 2 — irij — E — —jt for E 3> rrij . (6) 

The latter inequality indicates the relativistic limit of massive neutrinos which is the 
relevant limit for all neutrino oscillation experiments. 
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Since the state \u a , x) has flavour a at x = 0, we compute the probability to find the 
flavour p at x = L by \{vf3\v w L)\ 2 . With the explicit form (jSJ) of \v a , x) and with Eq. © 
we readily derive the standard formula of the neutrino transition or survival probability: 



P Va ^AL/E) 
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(7) 



It is easy to convince oneself that in the case of antineutrinos one simply has to make 
the replacement U — > U* in Eq. (J7J). 

Looking at Eq. (J7J), we read off the following properties of P Va ^ Up : 

• It describes a violation of family lepton numbers; 

• It is a function of mass-squared differences, e.g. of m 2 , — m\ and mj — m 2 ; 

• It has an oscillatory behaviour in L/E; 

• The relation P Va -> V/3 = Pp /3 ^p a is fulfilled as a reflection of CPT invariance; 

• P Ua -* v „ is invariant under the transformation U a j — > e l ^ a U a j e 1 ^ with arbitrary 
phases <p a , (fij. 

Specializing Eq. (|7J) to 2-neutrino oscillations with flavours a ^ /3, one obtains the 
particularly simple formulas 

U=( C ° S l Sin ^, (8) 
V — sin 6 cos v I 

P Ua ^ U0 = sin 2 29 X Ml - cos , (9) 

Pua^U/3 = Pvp—>V a ) (10) 

Pv a -*v a — Pvp^vp — 1 ~~ Pv a ^Vj3 • (11) 

In Eq. (JHJ) the rephasing invariance of P Va ^, v& has been used to reduce [/ to a rotation 
matrix. For 2-neutrino oscillations the probabilities for neutrinos and antineutrinos are 
the same. 

The oscillation phases in Eq. (J7|) have the form 



Am 2 L f Am 2 \ ZlMeV\ 

2.53 x 



2E \ leV / V E 



where Am 2 is one of the mass-squared differences. Alternatively, the oscillation phase is 
expressed by the oscillation length L osc through 

Am 2 L L ( E \ /leV 2 \ , . 

— S 2.— with L« c = 2.48,„ x (_) x (_). (13) 

In a neutrino oscillation experiment, the sensitivity to Am 2 is determined by the 
requirement that the phase (|12j) is of order one or not much smaller; this requirement 
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depends on the characteristic ratio Lj E relevant for the experimental setup. By definition, 
SBL experiments are characterized by a sensitivity to Am 2 > 0.1 eV 2 , while terrestrial 
experiments with a sensitivity below 0.1 eV 2 are called LBL experiments. Note that 
using longer baselines or smaller energies moves the sensitivity to smaller mass-squared 
differences (see Eq. (1T2J)). 

No neutrino oscillation experiment has found a mass-squared difference Am 2 ^0.1 
eV 2 , apart from the LSND experiment J7j. This experiment, however, has not not 
been confirmed by any other experiment and requires a forth neutrino which has to 
be sterile because the invisible Z width at LEP accommodates exactly the three active 
neutrinos. A sterile neutrino is defined as a neutral massive fermion which mixes with 
the active neutrinos, but has negligible couplings to the W and Z bosons. The result 
of the LSND experiment is at odds ^H] with the combined data sets of either all other 
SBL experiments or the solar and atmospheric neutrino experiments, when interpreted in 
terms of 4-neutrino oscillations. See, however, Ref. |19j . which shows that a 5-neutrino 
interpretation might reconcile the LSND result with the negative results of the other SBL 
experiments. 

Let us discuss some concrete examples for the sensitivity to Am 2 . For reactor neutrinos 
the energy is about E ~ 1 MeV. Therefore, with L ~ 1 km, the sensitivity is Am 2 ~ 10~ 3 
eV 2 ; this is the case of the CHOOZ [20] and Palo Verde [21] experiments. On the other 
hand, the KamLAND reactor experiment [7] with L ~ 100 km is sensitive down to 
Am 2 ~ 10~ 5 eV 2 . Experiments exploiting the source of atmospheric neutrinos with 
typical energies E ~ 1 GeV have as maximal baseline the diameter of the earth, i.e. 
L < 10 4 km, from where the sensitivity Am 2 > 10~ 4 eV 2 follows. Solar neutrinos have a 
very long baseline of L ~ 150 x 10 6 km and rather low energies of E ~ 1 -v- 10 MeV. Thus, 
in principle with solar neutrino experiments one can reach Am 2 ~ 10~ n eV 2 . 

Quantum-mechanical aspects of neutrino oscillations: Our phenomenological de- 
rivation of the oscillation probability ((Zj) needs some quantum-mechanical support (for 
extensive reviews see Refs. [221 [23J). The points of our derivation which need justification 
are the following: 

• We have employed stationary neutrino states, i.e., the plane waves have the same 
energy E but different momenta pj; 

• The usage of plane waves for the neutrino mass eigenstates raises the question how 
this is compatible with the localization of neutrino source and detection in space. 

In order to find limitations of the canonical oscillation probability we have to use a 
mathematical picture which is as close as possible to the actual situation of an experiment. 
Neutrinos are never directly observed; thus, for the description of neutrino oscillations it 
is required to consider the complete neutrino production-detection chain and to use only 
those quantities or particles which are really observed or manipulated in an experiment 
[21]. Consequently, we are lead to consider neutrino production and detection as a single 
big Feynman diagram where both the source and detector particle are described by wave 
packets which are localized in space, whereas the neutrinos are presented by inner lines 
of the big diagram; such models are called external wave packet models [23 123 and the 
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D particle 

Figure 1: The particles participating in the detector reaction at xd- 

big diagram is treated in a field-theoretical way [2H1- However, due to the macroscopic 
distance between source and detector the neutrinos are described as on-shell although 
they are on an inner line of the Feynman diagram 27 . 

We start with a consideration of the detector particle which is central to the first 
issue above [2H1 1213 EOj- In all experiments performed so far the detector particle is a 
bound state with energy Ed, either a nucleon in a nucleus or an electron in an atom. 
Such states are energy eigenstates. In addition, the detector bound state will experience 
thermal random movements with a whole spectrum of energies related to the temperature. 
However, there are no phase correlations between the different energy components of the 
thermal energy distribution which is, therefore, very well described by a density matrix 
which is diagonal in energy and the summation over this energy distribution is performed 
in the cross section corresponding to the big production-detection Feynman diagram 
|3U| I31j. In other words, in the amplitude corresponding to the big diagram, we have a 
definite detector particle energy E D + ^thermal; with an incoherent summation over the 
distribution of £^ ermal m the cross section. Concerning the n D final states in the detector 
reaction (see Fig. [TJ, we note that energy/momentum measurements are performed and 
that it is again summed incoherently over these energies E' Dk in the cross section. Thus, 
looking at Fig. ^ we come to the conclusion that in the amplitude the neutrino energy E 
is fixed and given by 

E v = Y1 E Dk ~ E D - dermal ■ 
k 

Since the oscillation probability (JJJ) is directly derived from the cross section, we come to 
following conclusion which has, in particular, been stressed by H. Lipkin [2T|l2*Hl 13 U| l3~T]: 



Neutrinos with the same energy E but different 
momenta pj = ^ E 2 — m 2 are coherent. 

Moreover, summation over the neutrino energies E is effectively performed in the cross 
section, i.e., it is an incoherent summation and no wave packets are associated with 
neutrinos of definite mass. 

Now we address the justification of the second point. Whereas it is necessary to assume 
localization in space of the neutrino source and detector wave functions, we stress that the 
final states in the source and detection reaction can be taken as plane waves; the reason is, 
as mentioned above, that the measurements performed with them are energy /momentum 
measurements and for calculating the cross section corresponding to the big Feynman dia- 
gram one has to sum incoherently, i.e. in the cross section, over the regions in phase space 
subject to kinematical restrictions according to the experiment. Denoting the detector 
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Figure 2: The wave function of the detector particle. The vertical line indicates the 
momentum difference |p 2 — Pi\- 

particle wave function in phase space by iPd(p), a reasonable range for the width of this 
function is given by ctd ~ 1CT 3 -j- 100 MeV, where these limits are given for a detector 
particle bound in an atom or in a nucleus, respectively In Fig. |2J we have symbolically 
drawn a detector particle wave function. Looking at this 1-dimensional picture it is easy 
to find the condition for the existence of oscillations: Assuming for simplicity only two 
neutrinos, the neutrino mass eigenstates have the same energy E but different momenta 
Pi, P2 ©• Suppose that we fix the momenta of the final states of the detector process 
by p' Dk (k — 1, . . . , n D , see Fig. |TJ). By momentum conservation, the values of the detec- 
tor wave function relevant for the amplitude of the complete production-detection cross 
section are given by t/jp (j^kPrn ~ Pj^) f° r eacn neutrino mass eigenstate, where £ is a 
unit vector pointing from the source to the detector particle. Let us assume now that 
J2kPDk ~ Pi^ — 0- 111 order to obtain coherence between the neutrino states with mass 
rrii and m 2 , the momentum p 2 must not fulfill \p 2 — Pi\ 3> <7d, because in that case we 
have (j2kPDk ~ P2^ — 0, the two mass eigenstates do not interfere and there are no 
neutrino oscillations. Therefore, we are lead to the condition 

\P2-Pi\ - Am 2 /2E<a D (15) 

for neutrino oscillations — see Fig. |2J Denoting the width of the detector wave function in 
coordinate space by <j x d, then with Heisenberg's uncertainty relation a x o ~ 1/2<7d the 
condition (JTHD is rewritten as |271 1231 132] 

(T xD < L osc /(An) . (16) 

This condition is trivially fulfilled because <j x d is microscopic whereas L osc is macroscopic! 
Similar considerations can be made for the neutrino source, with an additional condition 
taking into account that the source must be unstable (271 123 E3- Thus we can have 
confidence that the oscillation probability (|7J) holds for practical purposes. 

The present discussion is certainly not the most general one but, as we believe, it is 
reasonably close to practical applications. It does not include the discussion of a coherence 
length, i.e. a maximal distance between source and detector until which coherence between 
different mass eigenstates is maintained. There is a vast literature on that point, see for 
instance Refs. [2JH I2HJ OU EH] and citations therein. One detail can, however, be 
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immediately deduced from our discussion: The better one determines experimentally the 
energies of all final states in the detection reaction, the better one determines the neutrino 
energy E. Since the coherence length is approximately obtained by L osc E/AE [3^, where 
E is the average neutrino energy and AE the neutrino energy spread, it is evident that 
by pure detector manipulations one can influence the coherence length and, in principle, 
by infinitely precise measurements of the energies of the final particles in the detector 
process, one could make it arbitrarily long [2*^1 134j . 

Some steps in the above investigation of the validity of the oscillation probability (JJJ) 
have been criticized in the literature — see, for instance, Refs. (2HI EHl EHj — and even the 
external wave packet model is not undisputed [SHI- This suggests further investigations 
into the validity conditions of Eq. (|7|). However, since all attempts up to now to find 
limitations of (J7J) accessible to experiment were in vain, the standard formula for the 
probability of neutrino oscillations in vacuum seems to be very robust. 



2.2 Matter effects 

Matter effects fSl EE IHEj play a very important role in neutrino oscillations. Here we 
confine ourselves to the case of ordinary matter, which is non-relativistic, electrically 
neutral and without preferred spin orientation. It consists of electrons, protons and 
neutrons, i.e. fermions / = e~ , p, n. We denote the corresponding matter (number) 
densities by N e = N p and iV n . 

We want to present a simple heuristic and straightforward derivation of matter effects. 
First we notice that the vacuum oscillation probability (JJJ) can be derived by using the 
Hamiltonian 

H?* = ±UmW, (17) 

where rh is the diagonal matrix of neutrino masses. In order to obtain the effective 
Hamiltonian in matter, we have to add to ifj ac a term describing matter effects. With 
the above properties of ordinary matter we have the expectation values 

(hlM matter = ^N f 5^ . (18) 

Adding CC and NC contributions for neutrinos in the background of ordinary matter, the 
SM of electroweak interactions together with Eq. (|18|) provides us with the Hamiltonian 
density 

^ mat = % E "£(i-7sK 

V ^ a=e,fj,,r 

x E N f( 5 «f + T 3/l - 2 sin 2 9 w Q f ) , (19) 
/ 

where 5 a f comes from a Fierz transformation of the CC interaction, 1 $w is the weak 
mixing angle, and T%f L and Qf are the weak isospin and electric charge of the fermion /, 
respectively. Thus we arrive at the effective flavour Hamiltonian ^Sl QBj 

ttT = T^m^ + v^GVdiag (jV e - l -N n , ~N n) -^N n ) . (20) 



1 This term is non-zero only for v e interacting with the background electrons. 
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Then the evolution of a neutrino state undergoing oscillations and matter effects is cal- 
culated by solving 

do I 
i— = H^ at a with a= [ a M 1 . (21) 



dx 



where a e , a M and a T are the amplitudes of electron, muon and tau neutrino flavours. 

For active neutrinos {y a with a = e,fx, r) the neutron density N n can be dropped 
in the effective Hamiltonian, since it does not induce flavour transitions. Thus for active 
neutrinos only the electron density is relevant and there are no matter effects in 2-neutrino 

~~ * v t transitions. If sterile neutrinos exist, N n does affect v a —>■ v s because sterile 
neutrinos do not experience matter effects, i.e., N(u s ) = in H™ &t . 

If we consider 2-neutrino oscillations, we obtain from Eq. ()20j) the well-known 2-flavour 
Hamiltonian by subtraction of an irrelevant diagonal matrix and by using the 2x2 mixing 
matrix © [HZj: 

„t = J_ (A- Am 2 cos2# Am 2 sin 26 \ 
v AE\ Am 2 sin 26 -A + Am 2 cos 26 J K > 

with 

A = 2V2G F E [N(u a ) - N(u p )} (23) 

and N{u e ) = N e -N n /2, iV(i/ M ) = N(u T ) = -N n /2, N(y 8 ) = 0. This equation shows once 
more what we have discussed above. From now on we do not discuss sterile neutrinos 
anymore in these lecture notes. 

Performing all the analogous procedures for antineutrinos, we obtain if™ at by making 
the replacements U — > U*, matter term — > —matter term in if™ at of Eq. (|20p. 

In order to have an idea of the strength of matter effects, we estimate the "matter 
potential" A for the two most important cases, the sun and the earth: 



A = 2V2G F EN e 



core of the sun: 
1.5xlO- 5 eV 2 (jJfe), 
earth matter: 

2-3xlO- 7 eV 2 (^) (^f-) 



(24) 



The matter potential is the quantity which has to be compared with Am 2 in order to 
estimate the influence of matter on neutrino oscillations. 

In the evolution equation (|21j) the Hamiltonian is in general a function of the space 
coordinates through the electron density N e (x). Looking at the 2-neutrino Hamiltonian 
(jUJ), a MSW resonance [THj occurs if there is a coordinate x res such that 

A(x ICS ) = Am 2 cos 20 (25) 

holds. In this case, the probability of flavour transitions can be very big, in particular, if 
the neutrino goes adiabatically through the resonance point, as we shall see in the next 
subsection. 

In the following, we will adopt — without loss of generality — the convention Am 2 > 
and 0° < 6 < 90° in the Hamiltonian 122(1 . Therefore, the resonance condition is fulfilled 
for 6 < 45° and a suitable electron density. 
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i? = arctan(p/£) 



z 



Figure 3: Scattering at a slab of matter of thickness d of an incident wave propagating in 
z-direction 

On the validity of the effective flavour Hamiltonian: The heuristic derivation 
of matter effects resulting in the Hamiltonian (|2U|) did not allow to see under which 
conditions Eqs. (J20|) and (J2~T|) hold. We now want to discuss a method which gives us 
some insight into the limitations of these equations. We follow Refs. jHHl EHj • Suppose 
we have an incident wave propagating in ^-direction which traverses a slab of matter of 
"infinitesimal" thickness d. Then the incident wave undergoes multiple scattering at the 
scatterers / with number density Nf in the slab. According to Fig. |3] we sum over all 
scattered waves hitting a point on the z-axis at a distance I from the slab and obtain 



In this equation, A Ua j is the scattering amplitude. Eq. (j^B} can be considerably simplified. 
Firstly, it is reasonable to assume that in the term in the second line of this equation we can 
drop the limit p — ► oo because that contribution averages out. Secondly, for l/(Et) <C 1 
and a smooth behaviour of the scattering amplitude, we can also drop the term in the 
third line of Eq. ()26|) . Now we go one step further and consider a slab of finite thickness 
D = kd, where d is "infinitesimal" whereas k is a large integer. With the above-mentioned 
approximations we compute the phase of incident + transmitted wave through the finite 
slab as 



The expression in the curly brackets of the first line of this equation, where the incident 
wave is represented by 1 + iED/k ~ exp(iED/k), corresponds to the "infinitesimal" 
slab of thickness d = D/k and £ = (see Fig. EJ); the power k indicates that the wave 
propagates through k such "infinitesimal" slabs. This consideration yields the well-known 




(26) 




(27) 



exp {iED [1 + 2irN f A Va f(E,'& = 0)/E 2 }} . 
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V 




n 




2 sin 2 6 w + l/2 


—2 sin 2 6 1 ,,, 


+ 1/2 


-1/2 




2 sin 2 6 W - 1/2 


-2 sin 2 6 W 


+ 1/2 


-1/2 



Table 1: Vector coupling constants of the Hamiltonian density (J29|) . 



"index of refraction" 

n(u a ) = l + —J2 N f A Uaf (E, = 0), (28) 

where only the forward-scattering amplitude, i.e. i? = contributes. In n{v a ) we have 
indicated that different types of scatterers can be distributed in the slab. Furthermore, 
the derivation of the "index of refraction" is evidently not specific to neutrino scattering. 
The connection with neutrino scattering is made by the Hamiltonian density 

n Vaf = ^ p a7p (i - 75K FY (g { v' f) - g^-r*) / . (29) 

The coupling constants gy 1 '^ are found in Table H The axial- vector coupling constants 
do not contribute for ordinary matter (see Eq. 1)180) . As mentioned before, in the v e e~ 
scattering amplitude, there are CC + NC contributions, otherwise only NC interactions 
contribute. With the weak forward-scattering amplitude |3H1 

A„ af (E^ = 0) = ^-g^\ (30) 

Eq. (|2*7jl and Table UJ we compute 2ir X)/ NfA Ua f{E, i? = O)/^ and obtain exactly the 
matter terms as in 7i™ at in Eq. 1)20)1 . 

Now we use this second derivation of the matter potential as a means to check the 
validity of the Hamiltonian 1)271)1 |37J] . We define two lengths: rf scat t is the average distance 
of scatterers in ordinary matter, <i var ~ Nf(x)/\ dN Vy | is the typical distance of matter 
density variations. Then we note the following conditions for the validity of Eqs. (|2*0j) and 

* Ed scatt > 1 A^ Br °g lic < 27r4 



* 



'scatt 1 ' z,/lu scatt i 

*^scatt *^var j 
^scatt -^osc • 



The 2nd and 3rd condition arise from the requirement of describing the neutrino state 
evolution by the differential equation (|2~T|) . They are trivially fulfilled because d scat t is 
microscopic, whereas <i var and L osc are macroscopic. For instance for the sun we have 
(i V ar ~ 10 4 km. 

The first condition arises from the dropping of the term in the third line of Eq. ()26p. 
The quantity ^deBrogiie _ 2n / E is the de Broglie wave length of the neutrino. First we 
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have to estimate <i sca tt- It is easy the convince oneself that to a very good approximation 
the electron density in matter is given by 

N e ~Y e N A (p/lg), (31) 

where Y e is the number of electrons per nucleon, p is the matter density in units of g/cm 3 
and N A — 6.022 x 10 23 is Avogadro's constant. Let us check the first condition for the 
sun. According to the Solar Standard Model |Hj, in the solar core one has p ~ 150 and 
Y e ~ 2/3, from where it follows that N e ~ 100 N A cm" 3 and d scaU = N~ 1 ^ ~ 0.25 A. 
Thus the first condition above is reformulated as E d scat t — 130 £7/ (1 MeV) ^> 1. For solar 
neutrinos with E > 1 MeV this condition is fulfilled. For earth matter this condition holds 
as well, because c4 ca tt is a little larger than in the solar core. 

For other derivations of matter effects see Ref. j3H| and references therein. 



Survival probability for solar neutrinos: We want to conclude this subsection by 
some general considerations concerning solar neutrinos within the framework of 2-neutrino 
oscillations. A neutrino produced in the solar core traverses first solar matter, then travels 
through vacuum to the neutrino detector on earth; during the night, the neutrino has to 
traverse, in addition, some stretch of earth matter. We denote by P^ 2 the probability 
for |z/ e ) — > 1^1.2) in the sun. Furthermore, the probability of (^1,2) —> \v e ) transitions in 
earth matter is called Pf 2e - Note that with the 2x2 mixing matrix U of Eq. (JBJ) during 
the day we have 

pE = cos ^0 anc i pz = s m 2 9, (32) 

where 9 is the solar neutrino mixing angle. Then the survival probability of solar electron 
neutrinos is written as [4*T] 

p pS pE 1 pS pE I 

r "e-*Ve — r el r le r e2 r 2e ' 



2 y/PlPgP&Pi cos(5 m + Am 2 L vac /2E) . (33) 

In this formula, L vac ~ 150 x 10 6 km is the distance between the surface of the sun and 
the surface of the earth along the neutrino trajectory, i.e. Am 2 L vac /2E is the phase the 
neutrino acquires in vacuum. The phase acquired in matter is denoted by 8 m . 

Eq. (J33|) is completely general. One obvious questions arises: Under which conditions 
can the interference term in P Va _+ Ue be dropped? To investigate this point one can make a 
rough estimate. The size of the neutrino production region is approximately the diameter 
4ore °f the solar core. Therefore, for Am 2 £ corc /2E > 2ir with £ COTe ~ 10 5 km, we obtain 
the condition 

Am 2 /E > 10~ 14 eV (cos(5 m + Am 2 L vac /2E)) = . (34) 

This argument is not completely correct because of the matter effects, but numerical 
calculations give the same result |4*T] . On the other hand, experiments cannot measure 
the neutrino energy E with infinite precision and energy averaging in the vacuum phase 
occurs. If (Am 2 L/2E)(5E/ E) > 2n, where 5E is the uncertainty in the measurement of 
E, with an optimistic assumption of 5E/E ~ 10~ 2 we again arrive at Eq. (|34j) . 
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Figure 4: Allowed regions in the tan 2 6-Am 2 plane for the 2-neutrino oscillation solutions 
of the solar neutrino deficit. The upper region represents the LMA MSW solution, the 
lower one the LOW solution. This plot showing the status after the SNO results but 
before the KamLAND result has been taken from Ref. |42| . 



For the LMA MSW oscillation solution of the solar neutrino puzzle with Am 2 ~ 1CT 5 
eV 2 , condition (JMj) is very well fulfilled, and for the LOW solution — named after its "low" 
mass-squared difference of Am 2 ~ 1CT 7 eV 2 — this condition is still fulfilled 0]] . For the 
allowed regions in the tan 2 6-Am 2 plane of these solutions see Fig. 0] Note that this 
figure is relevant for the situation after the SNO results but before the KamLAND result 
which has ruled out the LOW region in the plot. For the recent history of solar neutrino 
oscillations see Ref. j3]. For the LMA MSW solution the resonance condition (j2HJ) plays 
an important role, as we will see in the next subsection. 



2.3 Adiabatic neutrino evolution in matter 

Now we want to discuss adiabatic neutrino evolution in matter. We focus on solar neutri- 
nos and confine ourselves to 2-neutrino flavours. Solutions of the 2-neutrino differential 
equation 

"Z) = " r ( Z ) (35) 

can be found by numerical integration and, for instance, the plots of allowed tan 2 6-Am 2 
regions of solar neutrino oscillations like Fig. H] are usually obtained via such numerical so- 
lutions, taking into account averaging over the solar neutrino production region. Eq. 
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refers to v e — *> v x transitions. In the next subsection, in the context of three neutrino 
flavours, we will discuss to which u^-u r flavour combination the amplitude a x refers. The 
Hamiltonian if™ at is given by Eq. (J2*2*|l . 

In general, solutions of Eq. (}35|) will be non-adiabatic. However, since the result of 
the KamLAND experiment we know that the solar neutrino puzzle is solved by neutrino 
oscillations corresponding to the LMW MSW solution; we will argue here that this solution 
behaves very well adiabatically. 

For the consideration of adiabaticity (see, for instance, Ref. |43j) we need the eigen- 
vectors of H™ at , which are defined Vx by 

H™\x)*l> mj {x) = Ejixtynjix) , (36) 

where 

/ cos6» m \ ( sm9 m \ , s 



In the 2-flavour case discussed here and with the real Hamiltonian (|22jl. one parameter, 
the matter angle 9 m , is sufficient to parameterize the two eigenvectors (p?7|) . With Eq. (J2*2*j) 
the matter angle is expressed as 

tan 20 

tan26 m (x) = jj-^r with A(x) = 2V2EN e (x) . (38) 

I _ 

Am 2 cos 29 

Note that for A — > 0, i.e. negligible matter effects, we obtain 9 m —>■ 6, i.e., the matter 
angle becomes identical with the (vacuum) mixing angle of U of Eq. (JHJ) • 
The full solution of the differential equation (|33|) can be written as 

iP( x ) = a j (x)^ mj (x)e~ ilp ^ x) with <pj = f* dx'E^x') . (39) 

j=l,2 Jx 

Having in mind solar neutrinos, we identify Am 2 with the solar mass-squared difference 
and use the initial condition that at Xq a v e is produced. This leads to 

*(*.) = ( i) =*■ ( ai< f°\ = C0S9 A X °\ ' (40) 

Eq. (}39j) is completely general. A solution is called adiabatic if the following condition 
is fulfilled. 

Adiabaticity : a\^{x) ~ constant. (41) 
In that case we have a 1]2 (x) — 01,2(^0) an d with Eq. (f4*Uj) we obtain 

Pei = cos 2 9 m (x ) , P| = sin 2 6 m (x ) . (42) 

Using now that the interference term in Eq. (|3*3*|) can be dropped for the LMA MSW 
solution and taking into account Eq. (J32J), we compute 



P Ue ~,v e = cos 2 9 m (x ) cos 2 6 + sin 2 9 m (x ) sin 2 6 

= -(l + cos2fl m (x )cos2#). (43) 
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This is the well-known survival probability for adiabatic neutrino evolution from matter 
to vacuum. Since we used Eq. (J3~2"j) . it does not include earth matter effects. 

Let us now consider effects of non-adiabaticity. Plugging the full solution into 
Eq. (JHEJ), we derive an equation for the coefficients ai j2 : 



d_ ( ax \ = ( -^e* ^ ( ai 
dx{ a 2 J { efce-* 



«2 



(44) 



where we have used the definition 

( p = (p 1 - ( p 2 = f dx'AE with AE = E x - E 2 . (45) 



With the Hamiltonian (1221) we calculate 



AE = — y(A — Am 2 cos 29) 2 + (Am 2 sin 26) 2 . (46) 

If along the neutrino trajectory the phase factor exp(i(p) in Eq. (J44)) oscillates very often 
while the change in 9 m is of order one, then di^ will be approximately zero. This suggest 
the introduction of the "adiabaticity parameter" 



7(a) = AE/2\9 m \ . (47) 

Therefore, we conclude 

Adiabaticity 7 > 1 (48) 

along the neutrino trajectory. In vacuum, 7 = 00, in agreement with a 12 being exactly 
constant. Non-adiabaticity is quantified by the probability P c for crossing from ip m i to 
It corrects Eq. flM)) to 45J 

P Va -+u e = \ + - Pc) cos 26 m (x ) cos 29 . (49) 

Again, as in Eq. (|4*H|) . this form of the survival probability holds where Eq. (pMjl is valid. 
The adiabaticity parameter 7 can be used to find a mathematically exact upper bound 
on P c 

Reverting to solar neutrinos, we use that the electron density in the sun fulfills 

N e (x) oc exp (— x/tq) , (50) 

except for the inner part of the core and toward the edge jH], with ro — R/ 10.54 ~ 
6.6 x 10 4 km <-> 3.3 x 10 20 MeV -1 ; R is the solar radius. We want to estimate 7 for the 
case that the neutrino goes through a resonance (|25jl in a region where the exponential 
form of the electron density provides a good approximation. Then adiabaticity will be 
rather well fulfilled if it is fulfilled at the resonance point |4^1 1431 14^ ITTj . There, 7 is given 
by 

Am 2 sin 2 29 AmVn 9 „ , . 

7res = , ~ x cos 29 tan 2 29 . 51 

2E cos 29(\ A\/A) res 2E v ; 



16 



For large 9 and close to the solar edge this estimate for adiabaticity is not correct (see 
Ref. @E]). 

Now we turn to the characteristics of the LMA MSW solar neutrino solution with 
Am 2 ~ 7 x KT 5 eV 2 and 9 ~ 34° [HE]. 2 Writing the oscillation length as 

, (1 x 10- 5 eV 2 \ / E \ 
L ~ = 35h °( Am 2 J(nfev)' (52) 

it is evident that in this case the flavour transition occurs inside the sun. Furthermore, 
with Eq. (|5T|) we estimate 7 res ~ 10 3 and the LMA MSW oscillation solution is clearly 
in the adiabatic regime. Other interesting features of the LMA MSW solution emerge by 
considering the limit of increasing neutrino energy: 



E T 



x res moves outward, 

P s ^ o ps^i V 6 ) 

Pv^v e -+±(i-cos20) = sin 2 0. 



The first line follows from Eq. (|25jl and the fact that N e (x) is monotonically decreasing 
from the center to the edge of the sun. 3 The other properties in Eq. (|53|) are read off from 
Eqs. (jHEJ), P^j) and (jUIJ), respectively. We remind the reader that the latter limit does not 
include earth matter effects, therefore, it holds only during the day, just as Eq. (J33)). Note 
that from P^ 2 -> 1 we conclude that for large enough solar neutrino energies we have a 
transition 

|^e)core ~~ > l^edge • (54) 

In Fig. 03 the solar v e survival probability is depicted. For historical reasons, P Ve -*u e 
is depicted not only for the LMA MSW, but also the LOW and vacuum oscillation 4 
probabilities which are both ruled out now. The dotted lines refer to daytime, i.e., when 
earth matter has no effect. The LMA MSW best fit of Ref. [IHj is given by Am 2 = 
5.0 x 10~ 5 eV 2 and tan 2 6> = 0.42 (after the SNO but before the KamLAND result). Let 
us now make a small numerical exercise. From the best fit we get sin 2 9 = 0.30, which 
should give P Ue -*u e for large E. Indeed, looking at the right end of the dotted curve in the 
right LMA panel where E — 15 MeV, we find excellent agreement. We see that neutrino 
energies around 10 MeV are already "large" in the sense of the limit (J53|) . On the other 
hand, let us take the limit E — > 0. In that limit the neutrino does not go through a 
resonance because Eq. ()25|) cannot be fulfilled. Therefore, because of the short oscillation 
length we expect an averaged survival probability P Ue ^ Ue = 1 — | sin 2 29 (see Eqs. JHJ) 
and (JUJ)). 5 Numerically, we obtain 1 — | sin 2 29 = 0.58, in excellent agreement with the 
left ends of the LMA panels. We have thus demonstrated here that qualitatively the LMA 
MSW solution can be quite easily understood. 



2 Before the result of the KamLAND experiment the mass-squared difference was a little lower at 
Am 2 - 5 x 1(T 5 eV 2 — see, e.g., Refs. [12 EE]. 

3 Note that a neutrino can go twice through a resonance if it is produced in that half of the sun which 
looks away from the earth. 

4 Here the oscillation length is of the order of the distance between sun and earth. 

5 We thank D.P. Roy for a discussion on this point. 
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Figure 5: Comparison of the survival probabilities for LMA MSW, LOW and vacuum 
solar neutrino oscillations for an electron neutrino created in the center of the sun. The 
figure has been taken from Ref. j3H]- The dotted lines refer to the daytime, whereas the 
dashed line includes v e regeneration during the night, when earth matter is effective. The 
full line refers to the average survival probability. 

We have not discussed earth matter effects here. Looking at the right LMA panel in 
Fig.0 we see that, at energies E^IO MeV, during the night when earth matter effects are 
operative the v e survival probability is a little larger than during the day. This introduces 
a small day-night asymmetry due to v e regeneration in earth matter. For a qualitative 
understanding of this effect see Ref. |5D] . 

Small Am 2 and the limit of strong non-adiabaticity: Here we want to discuss the 
question How do solar neutrinos approach vacuum oscillations in the limit of small Am 2 ? 
This is a non-trivial question because for small Am 2 there is always a point where the 
resonance condition (J25|) holds. Of course, this is a purely academic question since we 
know that solar neutrinos follow the LMA MSW solution for which the resonance does 
play a prominent role. It is, however, an interesting question from the point of view of 



18 



neutrino evolution in matter in general. 

When we speak of vacuum oscillations of solar neutrinos we mean mass-squared dif- 
ferences of the order of Am 2 ~ 10~ 10 eV 2 for which the oscillation length is of the order 
of the distance between sun and earth. From the resonance condition (}25|) we read off 
that in this case the resonance is very close to the solar edge (in this context, see N e (x) 
in Ref. [Sj). From Eq. (jHEJ) it follows that 9 m (xo) = n/2 and that 9 m = holds until 
very close to the resonance. It is easy to check that the width of the resonance, i.e. the 
distance where the change in A(x) is of the order of Am 2 cos 29, is roughly r , where the 
electron density N(x) behaves like Eq. (}50j) . Close to the edge of the sun the electron 
density drops steeper and the resonance width is smaller. Then one can check that while 
the neutrino crosses the resonance and also between the resonance and the solar edge one 
can approximate the phase ip of Eq. by a constant, which can be taken as its value at 
the solar edge. This is just the opposite of adiabaticity — see the discussion after Eq. (jUj) . 
For constant (p the solution of Eq. is given by 

ai(x)\ = ( cosA6> m -e ltfi sin A9 m \ ( ai(z )\ /rr\ 
a 2 (x)J \e-^smAe m cosA# m J \a 2 (x ) J 

with A6 m = 9 m {x) — 9 m (x ). In the case under discussion, the solution (|55j) is valid from 
just before the resonance till the solar edge and in the beginning of the vacuum. The 
initial conditions are 9 m (x ) = n/2 and ai(x ) = 0, a 2 (x ) = 1 . Then with Eq. at 
the solar edge, we get a% = e llf cos 9, a 2 = sin#. Plugging this result into the full solution 
(I39JI . we find what we expect: After the resonance, i.e. from the the solar edge onward, 
neutrinos perform ordinary vacuum oscillations; even the matter effects accumulated in 
the phases ip\ t 2 cancel. Inside the sun before the resonance, oscillations are completely 
suppressed because of the matter effects, due to A > Am 2 | cos 29 \. For further details 
see Ref. HE]. 

2.4 3-neutrino oscillations and the mixing matrix 

Neutrino mass spectra: Up to now all possible neutrino mass spectra, i.e. hierarchi- 
cal, inverted hierarchy, degenerate, etc., are compatible with present data. By convention 
we have stipulated mi < m 2 or Am\ x = Am 2 & > 0, where we have used the definition 
Am 2 k = m 2 — m\. Thus there are two possibilities for according to Am 2 x >0. Cor- 
respondingly, there are the two types of spectra depicted in Fig. EJ The "normal" and 
the "inverted" spectrum [21]. The hierarchical spectrum emerges in the limit m\ — ► of 
the normal spectrum, whereas the "inverted" hierarchy is obtained from the inverted 
spectrum with m 3 — > 0. 

In the 3-neutrino case the atmospheric mass-squared difference is not uniquely defined. 
If by convention we use for Am 2 tm the largest possible mass-squared difference, then for 
the normal spectrum we obtain Am 2 tm = Am 2 x and for the inverted case we have 
Am 2 tm = Am 2 3 . 

The mixing matrix: The most popular parameterization of the neutrino mixing matrix 
is given by [5^] 

U = U 23 U 13 U 12 (56) 
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Figure 6: The normal and the inverted 3-neutrino mass spectrum. 
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(57) 
(58) 
(59) 



The mixing angle 8u = # is probed in the solar neutrino experiments and in the Kam- 
LAND experiment. The mixing angle #23 = #atm is the relevant angle in atmospheric and 
LBL neutrino oscillation experiments. Up to now now effects of a non-zero angle #13 have 
not been seen. 

3-neutrino oscillations: The principle of 3-neutrino oscillations is a consequence of 
Amg <C Am^ tm and can be summarized in the following way. 



Aml tm L/(2E) ~ 1 AmlL/(2E) « 1 
Solar v oscillations frozen in atm. /LBL v osc. 



Am 2 e L/(2E) ~ 1 Aml tm L/(2E) » 1 
Atm. v oscillations averaged in solar v osc. 



Since, at present, #13 = is compatible with all available neutrino data, only an upper 
bound on this angle can be extracted. The two most important results in this context 
are: 



The non-observation of v e disappearance in the CHOOZ [zOj and Palo Verde 
experiments; 



The non-observation of atmospheric v e or v e disappearance [3]. 
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Also solar neutrino data have an effect on # 13 , although small, because in the 3-neutrino 
case the solar survival probability is given by 

^4€JH^i^3^) + 4, (so) 

where Pj%L, v is a 2-neutrino probability, calculated with the electron density cl 3 N e instead 
of N e . The result of fits to the neutrino data is 

sin 2 e 13 < 0.05 at 3a . (61) 

As evident from the parameterization 1)56)1 . CP violation disappears from neutrino oscil- 
lations for 6*i3 — > 0. 

In atmospheric neutrino experiments the mixing parameter which is extracted is 

sin 2 2^3~4|^ 3 | 2 (l-|^3| 2 ), (62) 

where we have taken into account the smallness of #13. The Super-Kamiokande experiment 
finds a result [S3] compatible with maximal mixing: 

sin 2 29 23 > 0.92 at 90% CL. (63) 

This leads to 0.60 < {U^l < 0.80. Note that 6 23 = 45° corresponds to \U^\ = l/y/2 ~ 
0.707. 

Now we want to discuss the states into which solar and atmospheric neutrinos are 
transformed. In the limit $13 — > or U\z — > 1 this discussion is most transparent. 
In that limit, in terms of neutrino states, atmospheric neutrino oscillations are given 
by |^ t ) — > \v T ), — > \v T ), while v e and z/ e do not oscillate because their oscillation 
amplitude is given by sin 2 2# 13 . In the context of Eq. ()35j) the question was raised to 
which neutrino state the amplitude a x belongs. For the evolution of neutrino states in 
the sun, the mass eigenstate |z/ 3 ) = S23 1 + C23 1 z^-r) is an approximate eigenstate of H™ &t 
because lAmJ-J A (see Eq. ()24jl ). Therefore, the initial state \v e )e must transform with 
probability -P®_,„ into the state orthogonal to \u 3 ), from where it follows that (without 
earth regeneration effects) 

K)o -> -C23W11) + S23K) • (64) 

Note that Eq. (|50jl has been derived by using the approximate eigenvector property of 
1 1/3), but corrections for non-zero #13 have been made. 

The primary goal of LBL neutrino oscillation experiments is the check of atmospheric 
neutrino results, observation of v T appearance and the oscillatory behaviour in L/E 
The demonstration of the latter property would give us the final proof that atmospheric 
Up and disappearance is really an effect of oscillations. Let us make a list what we 
would wish to gain from LBL (K2K with L = 250 km, MINOS and CERN-Gran Sasso 
with L = 730 km) and very LBL experiments: 

* Precision measurements of Am 2 tm and 623', 

* Measurement of 13 ; 
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* Verification of matter effects and distinction between normal and inverted mass 
spectra; 

* Effects of ArriQ and CP violation. 

We want to stress the difficulty of measuring CP violation in neutrino oscillations: Such 
an effect is only present if sin #13 7^ and if the experiment is sensitive to both mass- 
squared differences, atmospheric and solar, thus the measurements must be accurate to 
such a degree that they invalidate the principle of 3-neutrino oscillations mentioned in 
the beginning of this subsection. Furthermore, the earth matter background in (very) 
LBL experiments fakes effects of CP violation which have to be separated from true CP 
violation. 

Very LBL experiments are planned with super neutrino beams [22] and neutrino facto- 
ries [56J. Super neutrino beams arc conventional but very-high- intensity and low-energy 
(E ~ 1 GeV) beams, with the neutrino detector slightly off the beam axis; the JHF- 
Kamioka Neutrino Project with L = 295 km is scheduled to start in 2007 jSHI- Neutrino 
factories are muon storage rings with straight sections from where well-defined neutrino 
beams from muon decay emerge. According to the list above, measuring v e — > and 

~~ * transitions is of particular interest; note that #13 is the only angle in the neu- 
trino mixing matrix which has not been measured up to now. One can easily convince 
oneself that this task is tackled by measuring the so-called "wrong-sign" muons; e.g., 
/^storage ~~ * /^detector corresponds to v e — > v^. Neutrino factories are under investigation. Op- 
timization conditions for running such a machine are considered in the ranges 20 ^^^50 
GeV and L>3000 km. For further references see, e.g., jSZHSSHSHllSni and citations therein. 

3 Dirac versus Majorana neutrinos 

In this section we focus, in particular, on Majorana neutrinos. Apart from the vanishing 
electric charge, also most of the popular extensions of the SM suggest that neutrinos have 
Majorana nature. Majorana particles have several interesting features, which make them 
quite different from Dirac particles. Unfortunately, in practice, due to the smallness of 
the masses of the light neutrinos, it is quite difficult to distinguish between both natures: 
Neutrinoless /3/3-decay seems to be the only promising road so far. In any case, neutrino 
oscillations do not distinguish between Dirac and Majorana neutrinos, since there only 
the states with negative helicity enter, where this distinction is irrelevant. Though family 
lepton numbers must be violated for transitions v a — > z/g, the total lepton number remains 
conserved. 

3.1 Free fields 

With two independent chiral 4-spinor fields ipL,n one can construct the usual Lorentz- 
invariant bilinear for Dirac fields which is called mass term: 

Dirac: — m [^RipL + H.c.) = —m^ip with ip — ip L + ip R . (65) 
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With only one chiral 4-spinor ipi one obtains nevertheless a Lorentz-invariant bilinear 
with the help of the charge-conjugation matrix C: 

Majorana: -mi/i^C^ipi + H.c. = —-rmptp with ip — ipL + {iPl)°- (66) 

Note that from ipi one obtains a right-handed field with the charge-conjugation operation 
(i/jlY = C'JqiPl] in contrast to the Dirac case, this right-handed field is not independent 
of the left-handed field. 

The equation for free Dirac or Majorana fields in terms of the above defined spinors 
ip is the same: 

Dirac and Majorana: {i'y p d p — m)ip = . (67) 
In the formalism used here the Majorana nature is hidden in the 

Majorana condition: i\) = ip c . (68) 

Thus the solution of Eq. ()67j) is found in the same way for both fermion types. However, 
for Majorana neutrinos the condition ()68|) is imposed on the solution. This leads to the 
following observations. 

Dirac fermions: The field if) contains annihilation and creation operators a, b\ respec- 
tively, with independent operators a, b, and thus a Dirac fermion field has particles and 
antiparticles with positive and negative helicities. 

Majorana fermions: Because of the condition (J68|) the annihilation and creation operators 
are a, a\ respectively, and we have only particles with positive and negative helicities. 

The mass terms (|65J) and (JBTIj) are written in the way as they appear in the Lagrangian. 
The field equation (JBTj) is obtained from the Lagrangian by variation with respect to the 
independent fields. In this procedure the factor 1/2 in the Majorana mass term is canceled 
by the factor of 2 which occurs in the variation of the mass term because the fields to the 
left and to the right of C~ x are identical (see the first term in Eq. (J66j) ). 

Up to now we have discussed the case of one neutrino with mass m > 0. Let us now 
consider n neutrinos. Then in the Dirac case we have the following mass term: 

Dirac: — {vrM.vl + H.c.) = —vrhv . (69) 

Here, Ai is an arbitrary complex n x n matrix and the fields vl,r are vectors containing 
n 4-spinors. In order to arrive at the diagonal and positive matrix rh we use the following 
theorem concerning bidiagonalization in linear algebra. 

Theorem 1 If M. is an arbitrary complex nxn matrix, then there exist unitary matrices 
Ul,r with U r M.Ul = rh diagonal and positive. 

Applying this theorem we obtain the 

physical Dirac fields u' — v' L + v R with vl,r = Ul,r v'l r ■ (70) 

Note that the U(l) invariance of the mass term under u^r —>■ e %OL v^^R corresponds to total 
lepton number conservation, provided the rest of the Lagrangian respects this symmetry 
as well. 
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Switching to the Majorana case we have the following mass term: 

Majorana: ]^lC' l Mu L + H.c. = )^v' T L C~ l m v' L + H.c. = 

—-v'rhv'. (71) 
Now the mass matrix M. is a complex matrix which fulfills 

M T = M . (72) 

This follows from the anticommutation property of the fermionic fields and C T = — C. 
The diagonalization of the mass term proceeds now via a theorem of I. Schur |bT| . 

Theorem 2 If M. is a complex symmetric nxn matrix, then there exists a unitary matrix 
Ul with UJ^M-Ul = rh diagonal and positive. 

With this theorem we obtain the 

physical Majorana fields v' — v' L + (v'i) c with vl = Ulv' l . (73) 

The mass term (|71|) not only violates individual lepton family numbers just as the Dirac 
mass term (JBTJj) . but it also violates the total lepton number L = J2 a L a - 



3.2 Majorana neutrinos and CP 

It is easy to check that the Majorana mass term corresponding to the mass eigenfields, 
the second expression in Eq. (fTTj) (now we drop the primes on the fields i>l), is invariant 
under [TJl 

VjL Pj i Cv* L with p 2 - = 1 . (74) 

Thus Majorana neutrinos have imaginary CP parities. One can check that Vj = ujl + 
(z/jx) c transforms into pjiCv* under the transformation ()74|) . If we supplement this CP 
transformation by 

W+ - -e(l*)W- with e{p) = (-1)**+* , (75) 
U ^ ~Ct aL (a = e,/i,r) (76) 

and require that the CC interactions are CP-invariant then we are lead to 

U M = iU* M p with p= diag(p 1 ,p 2 ,Ps) ■ (77) 

We denote the mixing matrix for Majorana neutrinos by Um- From Eq. (J77J) we derive 

CP invariance U M = Oe^, (78) 

where O is a real orthogonal matrix. In the CP-conserving case the matrix O is relevant 
for neutrino oscillations since the phases in Um of Eq. (fTSj) drop out due to the rephasing 
invariance of the oscillation probability (see Section P2.1J1 . Note that one can extract an 
overall factor exp(z7r/4) or exp(— in/ 4) from the phase matrix in Um and absorb it into 
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the charged lepton fields. In this way, one obtains another commonly used form of Um in 
the case of CP conservation. 

For CP non-invariance one obtains 

U M = Ue l& . (79) 

The Majorana phases ctj cannot be removed by absorption into UjL because the Majorana 
mass term is not invariant under such a rephasing. However, one of the three phases can 
be absorbed into a the charged lepton fields, thus only two of the Majorana phases are 
physical. In neutrino oscillations only U ()56j) with the KM phase 5 is relevant. The two 
Majorana phases play an important role in the discussion of neutrinoless /3/3-decay and 
can be related to the phases which appear in leptogenesis and the baryon asymmetry 
from leptogenesis (521 (f° r reviews see, for instance, Refs. |64j). though in general the 
low-energy CP phases are independent of the phases relevant in leptogenesis [EHl EHj (see 
also Ref. [nij and citations therein). 



3.3 Neutrinoless /3/3-decay 

Up to now, the lepton number-conserving /3/3-decay (Z, A) —>■ (Z + 2, A) + 2e~ + 2v e has 
been observed in direct experiments with seven nuclides [68 . There is also an intensive 
search for the neutrinoless /3/3 or (/3/3)ou decay (Z, A) — > [Z + 2, A) + 2e~, where the total 
lepton number L = L e + + L T is violated (AL = 2). The distinct signal for such a 
decay is given by the two electrons each with energy E e = (M(Z, A) — MiZ + 2, A))/2, 
i.e. half of the mass difference between mother and daughter nuclide. In contrast to 
/3/3-decay with neutrinos, there is no unequivocal experimental demonstration for such a 
decay. The most stringent limit comes from 76 Ge with TfY 2 > 1.9 x 10 25 yr (90% CL) 169,. 
For reviews on ((3(3)q u decay experiments see Refs. (BE1EI]- 

There are several mechanisms on the quark level which induce two transitions n — > p 
in a nucleus without emission of neutrinos. Here we discuss some of the popular ones (see 
Ref. [7T] for a review of {(3(3)q v decay mechanisms). 

• Effective Majorana neutrino mass: This mechanism uses the Majorana nature 
of neutrinos and can schematically be depicted in the following way: 

d — » u + + u e 

d ^ u + e - + 1 y e } Wick contraction 

This Wick contraction can be performed because v c - = Uj and, therefore, the 
electron-neutrino field can be written as 

VeL = E^A/)^ = E( f/ A/)e J ^^C(^) T . (80) 
3 3 1 

Thus one actually contracts Vj with Uj Vj. In other words, the relevant neutrino 
propagator is given by the expression 

(0\Tis eL ( Xl )vJ L (x 2 )\0) = 

r d 4 n p-w-^-xi) 1 — 'Yk 

-£(^m) /V / ^- 2 — -^f^C. (81) 

V J (2n) 2 p 2 -m\ 2 
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One can show that it is allowed to neglect m| in the denominator of the Uj propa- 
gators [70], which leads then to the effective Majorana neutrino mass 

(m>=£(tfAr)i"V. (82) 
j 

For this mechanism of ((3(3)q v decay, the decay amplitude is proportional to \{m)\. 

The present lower bounds on Ty 2 f° r 76 Ge lead to an upper on | (m) | varying between 
0.33 and 1.35 eV, depending on the models used to calculate the nuclear matrix 
element [611] 172] . Future experiments want to probe down to a few x 0.01 eV 

M- 

In the case of CP conservation (see Eq. ([780) the effective Majorana neutrino mass 
is given by (m) = Oej Pj m j- For opposite signs of pj = ±1, i.e. opposite CP 
parities, cancellations in (m) naturally occur. 

Intermediate doubly charged scalar: In left-right symmetric models, doubly 
charged scalars occur in the scalar gauge triplets and induce ((3(3)q v decay via the 
mechanism [7T] depicted symbolically in the following way: 

d — > u + S~ \ 
d^u + S- / 



SUSY with R-parity violation: Within this framework, several mechanisms for 
inducing ((3(3)q v decay exist [H]. The gluino (g) exchange mechanism uses the 
Majorana nature of gauginos [73] : 



d 



d 



e + u 



u 



u + g 
u + g 



Wick contraction 



In this picture, u is an up squark. 

For the check of the Majorana nature of neutrinos, (/3/3)ou decay is most realistic possi- 
bility. Another line which is pursued is the search for v e from the sun [73] . Both types of 
experiments aim at finding |AL| = 2 processes. 

In view of several mechanisms for ((3/3) o v decay where some even do not require neu- 
trinos, the question arises if an experimental confirmation of this decay really signals 
Majorana nature of neutrinos. The question was answered affirmatively by Schechter and 
Valle [7^], and Takasugi [7B]. Their statement is the following: If neutrinoless f3(3-decay 
exists, then neutrinos have Majorana nature, irrespective of the mechanism for /3/3-decay. 
We follow the argument of Takasugi and show the following: If ((3(3)q v decay exists, then 
a Majorana neutrino mass term cannot be forbidden by a symmetry. 
Proof: We assume for simplicity that there is no neutrino mixing. As will be seen the 
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arguments can easily be extended to the mixing case. We introduce phase factors 77 and 
assume the symmetry 



e L -> 7] e e L , v eL -> rj u v eL , 

Invariance of the weak interactions under this symmetry give the following relations: 

VeLl P e L W+ rf v r) e r) W = 1 , (84) 
u L Yd L W+ => r]* u r] d r] W = 1 . (85) 

Now we assume the existence of (fl/3) 0u decay, which introduces a further relation: 

d L + d L ^ u L + u L + e L + e L => (VdVuVe) 2 = 1 (86) 

Eqs. (|H1J) and (J85J) together lead to r) v = rj d rj u rj e . Thus with Eq. we arrive at r)l — 1, 
and the mass term oc v1 L C~ x v eL cannot be forbidden by v eL — > rj u v eL . Q.E.D. 

The above argument does not necessarily show that a Majorana neutrino mass term 
must appear, but in general this will be the case for renormalizeability reasons: All terms 
of dimension 4 or less which are compatible with the symmetries of the theory have to be 
included in the Lagrangian. 



3.4 Absolute neutrino masses 

With neutrino oscillations, the lightest neutrino mass which we denote by m cannot be 
determined. Nowadays three sources of information on m are used: 3 H decay, {(3f3)o u 
decay and cosmology [77] . For an extensive review see Ref. |51|. 



Tritium and neutrinoless /3/3-decays: In the decay 3 H — > 3 He + e~ + u e , with an 
energy release Eq — 18.6 keV, the region near the endpoint of the electron recoil energy 
spectrum is investigated. Deviations from the ordinary (3 spectrum with a massless v e 
are searched for. Assuming that the v e has a mass m^, the Mainz and Troitsk experi- 
ments have both obtained an upper limit mp < 2.2 eV at 95% CL. The future KATRIN 
experiment plans to have a sensitivity mp > 0.35 eV jlj. 

Since we know that the electron neutrino field is a linear combination of neutrino mass 
eigenfields because of neutrino mixing (see Eq. ©), the 3 H decay rate is actually a sum 
over the decay rates of the massive neutrinos multiplied by |t/ e j| 2 - Thus the question 
arises for meaning of mp which is extracted from the data. Farzan and Smirnov [7H] 
have argued in the following way. Firstly, the bulk of data from which the bound on 
mp is derived does not come from the very end of recoil spectrum. Secondly, the energy 
resolution in 3 H experiments, even for the KATRIN experiment, is significantly coarser 

than \J Am^tm • Therefore, the three masses rrij can effectively be replaced by one mass 
mp and the different expressions for mp found in the literature, e.g. mp = J2j \U e j\ 2 rrij, 
are all equivalent. The simplest choice is mp = m Q (for other papers on this subject see 
citations in |78j). Thus the results of the 3 H decay experiments give information on the 
lightest neutrino mass m . 
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In order to make contact between m and ((3f3) 0v decay, the following assumptions are 
made: Clearly, one has to assume Majorana neutrino nature, but also that the dominant 
mechanism for {(30)q v decay proceeds via the effective Majorana neutrino mass 
(182)1 . which is then expressed as a function of mo and the parameters of the Majorana 
neutrino mixing matrix Um (JZHJ)- Note that this mechanism must be present for Majorana 
neutrinos, it might only be that it is not the dominant one. For the normal spectrum we 
have m = m\ and m = m 3 for the inverted spectrum (see Fig. |BJ). 



Normal spectrum: 

I HI 

Inverted spectrum: 

I HI 



m c{ 2 



ml + Am% s\ 2 e 



ah 



cf 3 + 



ml + AmL s 13 e 



(87) 



ml + Aml tm 



Am® c" , ( 



.2 

"12 ' 



+ 



m 



atm °12 ' 



c 2 13 + m Si 3 



In these expressions, (3\ and fa are general CP-violating phases which arise according to 

6 



the form of Um 
e V and J Am| 



For numerical estimates it is useful to have in mind that \ Am 



1 at in 



0.05 



0.008 eV. For references see, for instance, (ZHIIHIII and citations therein. 
In order to plot | (m) \ against m , the values of neutrino oscillation parameters are used 
as input. The phases fii^ are free parameters; for CP conservation the possible values of 
(3i y 2 are zero or tt. The plot of Fig. [7|has been taken from Ref. [7^]. In the upper panel, 
best fit values of the neutrino oscillation parameters are used, whereas in the lower two 
panels the uncertainty in these parameters is taken into account. The lower band on the 
left sides of these panels correspond to the normal spectrum, whereas the upper band to 
the inverted spectrum. 

What does this figure teach us? One can use experimental upper bounds on | (m) \ and 
m and compare with the allowed regions in the plot. It has been put forward in Ref. [EI] 
that such a comparison, even if {(3(3)q v decay is not found experimentally, might reveal the 
neutrino nature: If the minimum of \(m) \ predicted from mp and the oscillation parame- 
ters exceeds the experimental upper bound on |(77i)(, then the neutrino must have Dirac 
nature. There are, however, several difficulties with such a procedure: Other mechanisms 
for (/5/?)ot/ decay could destroy the validity of such a comparison; m$ will not reach the 
non-degenerate neutrino mass region even with the KATRIN experiment, and the same 
difficulty applies to the region with |(m)| below 0.01 eV; moreover, imprecise determina- 
tion of the oscillation parameters blur the picture as seen from the lower two panels of 
Fig. [7| However, the chances for this approach would be good, if the KATRIN experiment 
would find a non-zero m$. On the other hand, assuming the validity of the Majorana 
hypothesis, Fig. might give us information about the type of neutrino mass spectrum if 
|(m)| is measured or a sufficiently stringent experimental bound on it is derived |79j . 



3 They are simple functions of the Majorana phases and of 6. 
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Figure 7: Effective Majorana neutrino mass |(j7i)| versus lightest mass m . The figure 
is taken from Ref. jTD]. In the upper panel best fit values for the neutrino oscillation 
parameters have been used and only the phases fa and fa have been varied. (A non-zero 
best fit value for S13 is used in this plot.) In the lower two panels the uncertainty in 
the determination of the oscillation parameters has been taken into account. For the 
details and for the precise meaning of the lines in the shaded regions, which correspond 
to CP-conserving choices of the phases fa$, see Ref. [79*] . 

Neutrino masses and cosmology: A very interesting bound on the sum over all 
neutrino masses is provided by the large-scale structure of the universe |BZI and the tem- 
perature fluctuations of the cosmic microwave background (CMB) [83] . Usually, energy 
densities pi in the universe are given as fractions flj = pij p cr of the critical energy density 
p cr of the universe. Today's energy density of non-relativistic neutrinos and antineutrinos 
is given by (see, for instance, [84J) 

tt u h 2 = ^77^/(93.5 eV), (89) 
j 

where h ~ 0.7 is the Hubble constant in units of 100 kms -1 Mpc -1 . Of the three active 
neutrinos, at least two of them are non-relativistic today since J Am@ ~ 0.008 eV and 

\J Am^tm ^> \J Am®, while the neutrino temperature today corresponds to 1.7 x 10~ 4 eV 
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Roughly speaking, an upper bound on Vt v is obtained in the following way. Hot dark 
matter tends to erase primordial density fluctuations. The theory of structure formation 
together with data on large-scale structures (distribution of galaxies and galaxy clusters), 
mainly from the 2 degree field Galaxy Redshift Survey [82J, gives an upper bound on 
Qv/Qm, where Q m is the total matter density. On the other hand, data on the temper- 
ature fluctuations of the CMB can be evaluated with the so-called ACDM model (a flat 
universe with cold dark matter and a cosmological constant) and allows in this way to de- 
termine a host of quantities like h, Q m , (the baryon density), etc. Since the quantities 
extracted in this way agree rather well with determinations based on different methods 
and assumptions |SS1, the ACDM model is emerging as the standard model of cosmol- 
ogy jHSj- Evaluation of the recent CMB results of the Wilkinson Microwave Anisotropy 
Probe (WMAP) and of the large-scale structure data give the impressive bound [EB] 
fl u h 2 < 0.0076 (95% CL); with Eq. (|H9*|) this bound is reformulated in terms of neutrino 
masses: 

Y.mj < 0.7 eV (95%CL). (90) 

i 

For a single neutrino we have thus rrij < 0.23 eV, since neutrinos with masses in the range 
of a few tenths eV have to be degenerate. 

Note that the ACDM model gives the result Q m ~ 0.3 and Q\ ~ 0.7; the latter 
quantity is the energy density associated with the cosmological constant A. Furthermore, 
one finds £l u ?S 0.015 and Qb ~ 0.04 -=-0.05, thus the main contribution to Q m must consist 
of hitherto undetermined dark matter components. 

3.5 Neutrino electromagnetic moments 

The effective Hamiltonians for neutrinos with magnetic moments and electric dipole mo- 
ments (electromagnetic moments) are given by the following expressions: 

Dirac: H° m = h R \ ^ v v L F^ + H.c, (91) 
Majorana: Hf m = --i/JC _1 A o^vlF^ + H.c. (92) 

In the Majorana case, from (C _1 o" M!y ) T = C~ 1 <j^ v and the anticommutation property of 
i>l, it follows that A T = —A and only transition moments are allowed. In the Dirac case, 
the electromagnetic moment matrix A is an arbitrary 3x3 matrix. Usually, a decom- 
position of A is made into A = fi — id with fp = /i being the MM matrix and d) = d 
the EDM matrix. However, since the neutrino mass eigenstates are experimentally not 
accessible, the distinction MM/EDM is completely unphysical in terrestrial experiments. 
This holds also largely for solar neutrinos, though in the evolution of the solar neutrino 
state the squares of the neutrino masses enter, and one can shift phases from the mixing 
matrix to the electromagnetic moment matrix and vice versa, which makes the distinc- 
tion MM/EDM phase-convention-dependent. For a thorough discussion of this point see 
Ref. [El. 

Bounds on A from solar and reactor neutrinos are obtained via elastic ve~ scattering 
where, due to the helicity flip in the Hamiltonians (j9Tf and (J92|) . the cross section is the 
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sum of the weak and electromagnetic cross sections [8*7] : 



da da w da em da em a 2 n (I 1 \ 2 

dT = lT + ^F Wlth ^ = ^|It-eJ^- (93) 

In the most general form, the effective MM is given by [HE] 

/Xg ff = oLaUgl + AA f a + . (94) 

In Eq. (I93j) . T is the kinetic energy of the recoil electron, //g is the Bohr magneton, and the 
flavour 3-vectors a± describe the neutrino helicity states at the detector. The expression 
(|94|) is basis-independent, thus it does not matter in which basis, flavour or mass basis, 
the quantities A, a± and the Hamiltonians ("9*""j) and ("9*!"]) are perceived 86J. 

In the following we will consider reactor and solar neutrinos. If the detector is close 
to the reactor, one simply has (in the flavour basis) a_ = 0, a + = (1, 0, 0) T , as the reactor 
emits antineutrinos with positive helicity. For solar neutrinos, the effective MM will, 
in general, depend on the neutrino energy E. Bounds on the neutrino electromagnetic 
moments have been obtained from both reactor jHH] and solar neutrinos jH3 EH] • 

From now on, we concentrate on Majorana neutrinos, where the antisymmetric matrix 
A contains only three complex parameters (in the Dirac case there are nine parameters). 
Thus, A can be written as 

A Q/ 3 = £ Q /3 7 A 7 , (95) 

in either basis, flavour or mass basis. For the approximations used in the calculation of 
the effective MM for solar neutrinos and the LMA MSW solution, //lmaj we re f er the 
reader to Ref. [HI]. The result is given by 

^ma(^) = |A| 2 - IA2I 2 + P?i (|A 2 | 2 - \Ai\ 2 ) . (96) 

The probability is defined before Eq. (|32j) . For reactor neutrinos the effective MM is 
immediately obtained as 

a 2 IX I 2 i I \ I 2 — I A I 2 _i_ I A I 2 

A^reactor l^e/xl ~ l^erl MVI ' I "n 

= |A| 2 - c^lAil 2 - sf 2 |A 2 | 2 - 2s 12 c 12 |A 1 ||A 2 | cos 5' , (97) 

where the second line is derived from the first line by using A = ?7mA. The phase 5' is 
composed of 5' = arg(A^A 2 ) + oi<i — a± with the Majorana phases atj defined in Eq. (JTUJ). 
Eq. (|9*3j) defines the vectors A = (A a ) and A = (Aj) in the flavour and mass basis, 
respectively. The length of A occurs in both effective moments, Eqs. (j^Ujl and (|97|). Note, 
however, that 

|A| 2 = ^Tr(AU) |A| = |A|, (98) 

i.e., the length of the vector of the electromagnetic moments is basis- independent. 

In Ref. [HI] a bound on | A| is derived, i.e., all three transition moments are bounded by 
using as input the rates of the solar neutrino experiments, the shape of Super-Kamiokande 
recoil electron energy spectrum and the results of reactor experiments at Bugey (MUNU) 
and Rovno |S8| . As statistical procedure a Bayesian method is used with flat priors and 
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minimization of x 2 with respect to 6*12, Am , |Ai|, |A 2 |, 5', in order to extract a probability 
distribution for |A|. The details are given in Ref. j^I], where the 90% CL bounds 

f 6.3 x 10~ 10 yUB (solar data) , , . 

' ' \ 2.0 x 10~ 10 fi B (solar + reactor data) ^ ' 

are presented. 

Some comments are at order |91j : 

> We want to stress once more that the bounds in Eq. (|99j) apply to |A| 2 = |A 



12 



v e/j I 



|A MT | 2 + |A re | 2 = J2j<k \^jk\ 2 ■ Therefore, all transition moments are bounded in a 
basis-independent way. 

> In Eq. ([96)1. in the limit P^ — > 0, the quantity |A 2 | drops out of the effective MM. 
Therefore, for small P^, solar neutrino data become less stringent for |A 2 | and thus 
for |A| as well. This is the case for Super-Kamiokande data (see the first bound in 
Eq. ([99))) where P^ is small due to the relatively high neutrino energy (see Eq. ([53)1 ). 

> Reactor data give a good bound on |A 2 | and are, therefore, complementary to present 
solar neutrino data. 

> The BOREXINO experiment [H2] could improve the second bound in Eq. (|99|1 by 
nearly one order of magnitude because it is sensitive to relatively low neutrino 
energies where P^ ~ 0.5. 

4 Models for neutrino masses and mixing 
4.1 Introduction and scope 

For a start, we present the problems of model building for neutrino masses and mixing 
by asking appropriate questions and formulating answers, if available. 

1. Can neutrino masses and mixing be accommodated in a model? 

This is no problem; the simplest possibility is given by the SM + 3 right-handed 
neutrino singlets Ur + L conservation which allows to accommodate arbitrary masses 
and mixing of Dirac neutrinos in complete analogy to the quark sector. 

2. Why are neutrino masses much smaller than the charged lepton masses? 
There are two popular proposals for a solution to this problem: 

• Radiative neutrino masses; 

• Seesaw mechanism 



Can one reproduce the special features of v masses and mixing? 
Let us list the features one would like to explain: 

Fl #0 ~ 34° (large but non-maximal); 
F2 atm ~ 45°; 
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ew/TeV scale 


RGeq 


GUT scale 







gauge symmetry 




±SUSY 


no/abelian/non-abelian 
flavour symmetry 



Figure 8: Symbolical presentation of the variety of possibilities for model building. The 
upper two blocks show the two energy scales at which the symmetries, indicated in the 
three blocks below, could apply. RG eq stands for renormalization group equation, ± for 
with/without. 

F3 \U e3 \ 2 = s? 3 <0.05; 
F4 Am|/Am 2 atm ~ 0.03. 

We have used the obvious notation where the subscripts and atm refer to solar 
and atmospheric neutrinos, respectively. 

To explain the listed features is the most difficult and largely unsolved task. There are 
myriads of textures or models — for reviews see Ref. j^H]- One of the problems is that there 
are still not enough clues where to start model building. Symbolically, the difficulties of 
model building are presented in Fig. |HJ If the symmetries apply at the GUT scale, the 
renormalization group equation transports relations from the GUT scale to the low scale; 
this transport could contribute to generate some of the desired features — for a review see 
Ref. (HH]. It is completely unknown if the explanation of the features Fl-4 is independent 
of general fermion mass problem or not. In the following we will assume independence. 

The scope of the following sections is to discuss some simple extensions of the SM by 
addition of scalar multiplets and of right-handed neutrino singlets vr. All these extensions 
will yield Majorana neutrino masses. We will only discuss the lepton sector. 

4.2 The Standard Model with additional scalar multiplets 

In the SM there are only two types of lepton multiplets: 

SU(2) x 17(1) 
Dl \ Y = — 1 left-handed doublets 

£r Y = — 2 right-handed singlets 

Here Y is the hypercharge and the underlined numbers indicate the weak isospin of the 
multiplet. 
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By forming all possible leptonic bilinears one obtains the possible scalar gauge multi- 
plets which can couple to fermions [§Tj : 



R 



D L ®D L 





1 

2 



i 







doublet 
singlet 
triplet 
singlet 



The hypercharge in this table refers to the leptonic bilinear. In the following we will 
discuss extensions of the SM with these scalar multiplets (except the trivial extension 
where only Higgs doublets are added). 



The Zee model: This model 
has the Lagrangian 



is defined as SM with 20 + rj + and, therefore, 



C 



■ ■ + 



f aP D T aL C-HT 2 D pL r/+ - /i 0l0 2 77 + + H.c 



(100) 



where the dots indicate the SM part and terms of the Higgs potential which are not 
interesting for the following discussion. Note the antisymmetry of the coupling matrix: 

fa/3 = —f/3a- 

Since no fermionic multiplets have been added to the SM, the ensuing neutrino masses 
can only be of Majorana type. Thus the total lepton number L must be broken, otherwise 
neutrino masses will be strictly forbidden. Let us assign a lepton number to 77+ from its 
Yukawa couplings in Eq. (jlOOj) . Then we have the following list of lepton numbers of the 
multiplets of the Zee model: 



L 



D, 



h,2 V~ 



1 







(101) 



Thus, L is indeed explicitly violated by the term 0{0 2 ?7+ m the Lagrangian. Such a term 
can only be formed if two Higgs doublets are present (if 0i = 02 = 0, then <f><<j> = 0). In 
the Zee model a neutrino mass matrix M. v appears at the 1-loop level. 

A lot of work has been done on the restricted Zee model jlOOj where only one Higgs 
doublet, say 0i, couples to the leptons. In this particularly simple version one has M. v oc 
(( m a — mfyfap), where m a ^ denote the charged lepton masses. Note that without loss 
of generality we have chosen a basis where the charged lepton mass matrix is diagonal. 
The restricted Zee model is practically ruled out now because it allows only maximal 



I A 



solar mixing |101j : furthermore, it requires serious fine-tuning, namely m^|/ er | ~ rn; lJetlU 
finally, the smallness of the neutrino masses has to be achieved by \f a p\ ^10~ 4 . If both 12 
couple to the leptons, then non- maximal solar mixing is admitted [102J, the fine-tuning 
problem is somewhat alleviated, but the third point remains. For further recent literature 
on the Zee model see, e.g. Ref. j!03j 



The Zee-Babu model: 

the Lagrangian 

C = 



This model is defined as jHH I104j SM + 77+ + k ++ and has 



+ 



f aP Dl L C-Hr 2 D^ 



+h 



aP l aR 



c- 



k ++ — jirj i] k ++ + H.c. 



(102) 
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with a symmetric coupling matrix h a p = hp a . In addition to the assignments (jlUlj) . we 
have L{k ++ ) = —2. Again, L must be explicitly broken, which is now achieved by the 
//-term in Eq. (jl02|) . 

Since in this model there is only one Higgs doublet, L is still conserved at the 1-loop 
level (see previous section on the Zee model) and the neutrino mass matrix appears at the 
2-loop level: M. v oc frhih*rhef with / = (/a/3), h = (h a a) and = diag (m e , m M , m T ). 
This model has the following properties |lU5j . Since / is antisymmetric, the lightest 
neutrino mass is zero; the solar LMA MSW solution and, e.g., a hierarchical neutrino 
mass spectrum require the fine-tuning \h^\ : \h^ T \ : \h TT \ ~ 1 : {m^/m T ) : (m^/mj 2 ; all 
scalar masses are in the TeV range; in order to reproduce the neutrino masses inferred 
from atmospheric and solar data one needs \f a p\, \h a p\ ~ 0.1, thus neutrino masses are 
naturally small as a consequence of the 2-loop mechanism; finally, rare decays like r — ► 3/i, 
// — ► e7 should be within reach of forthcoming experiments. 

We note that from the Zee model and Zee-Babu model we have seen that models 
with radiative neutrino mass generation are prone to excessive fine-tuning because the 
hierarchy of the charged lepton masses works against the features needed in the neutrino 
sector. 



The triplet model: This model is defined by adding a scalar triplet A to the SM and 

has the Lagrangian 



C = ■■■ + 



\g a(3 D T aL C-Hr 2 ^D pL + H.c. 



-|M 2 TrAtA - (//0tA0 + H.c." 
The electric charge eigenfields and the VEV of the triplet are given as 



(103) 



respectively. Note that g a p = gp a . As in the previous two models, we can make the 
assignment L(A) = —2, then L is explicitly broken by the //-term in Eq. (|lU3j) . The 
original model without the //-term and spontaneous L breaking [106J is ruled out by the 
non-discovery of the Goldstone boson and a light scalar at LEP. This model leads to 
the tree-level neutrino mass matrix M. v = VT{g a /3)- The LEP data require \vt/v\ ^ 0.03 
|107j . where v is the VEV of the SM Higgs doublet. However, if the coupling constants 
g a p are about 0.01 -v- 0.1, then in order to obtain small neutrino masses the triplet VEV 
must be much smaller, namely Vt ~ 0.1 -j- 1 eV. There are two ways to get a small vt'- 
Firstly, one can assume that M, |/t| 3> v, then \vt\ — \n\v 2 /M 2 (scalar or type II seesaw 
mechanism) I109j : secondly, with M ~ v, |/t| C « one has \vj-\ ~ \^\ |108j . In some 
sense the mechanism for obtaining small neutrino masses in the triplet model is analogous 
to the seesaw mechanism [UHl EH (see subsequent subsection) since in both cases one has 
to introduce a second scale much larger (smaller) than the electroweak scale in order to 
generate small neutrino masses. 
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4.3 The seesaw mechanism 



The seesaw mechanism [931 1941 is implemented in the simplest way in the SM + 3 Ur + L 
violation. This is primarily an extension in the fermion sector of the SM. Note that one 
could also choose two or more than three right-handed singlets vr. The starting point is 
the Lagrangian 



C 



Ej (<no?\, + %$A) d l + h.c. 

+ (iz/gC- 1 M^, + H.c), 



(105) 



where the mass matrix Mr of the right-handed singlets must be symmetric due to their 
fermionic anticommutation property. The number n# of Higgs doublets is irrelevant for 
the seesaw mechanism. Defining the mass matrices 

M t=^H v Fi> M D = i=5>A J; (106) 

where Mi is the mass of the charged leptons and Mp is the so-called "Dirac mass matrix" 
for the neutrinos, the total Majorana mass matrix for all six left-handed fields up is given 
by PHI] 



The VEVs vj fulfill v = ^j\ Vl \ 2 + . . . \v nH \ 2 = (V^Gp)- 1 / 2 ~ 246 GeV. 

The basic assumption for implementing the seesaw mechanism is < tur, where 
m D,R are the scales of Mp t R, respectively. With this assumption, one obtains the mass 
matrix of the light neutrinos 

M u = -MlM^M D , (108) 

which is valid up to corrections of order (mp/rriR) 2 . The mass matrix of heavy neutrinos 
is given by Ml cwy = Mr. Diagonal izing the mass matrices by 

{U R ) ] M f Ui = m e , V T M U V = m , (109) 

the neutrino mixing matrix for the light neutrinos is given by 

u M = (u e L yv. (no) 

As discussed earlier, phases multiplying Um from the left are unphysical because they can 
be absorbed into the charged lepton fields. 

The seesaw mechanism contains three sources for neutrino mixing: Mg, Mp, and Mr. 
Therefore, it is a rich playground for model building. It can also be combined with 
radiative neutrino mass generation — for an example see next subsection. Let us discuss 
the order of magnitude of the scale itir. If we choose as a typical neutrino mass m v ~ 



A m atm ~ 0.05 eV and assume vrtp ~ m^ T , we obtain vrtR ~ 10 s -r- 10 11 GeV. On the 
other hand, if mp is of the order of the electroweak scale, then tur ~ 10 15 GeV and there 
could be a possibility to identify it with the GUT scale. 
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4.4 Combining the seesaw mechanism with radiative mass gen- 
eration 

Let us now consider the SM with two Higgs doublets 4>i,2, add a single field vr and allow 
for L violation I112j . We do not employ any flavour symmetry. 

In this case the Yukawa coupling matrices Ai^ and Mo are 1x3 matrices (vectors!). 
Here we give a precise meaning to the scale mo and identify it with the length of the 
vector Mo- L violation is induced by the Majorana mass term with mass Mr > of 
the neutrino singlet vr. Then at the tree level the seesaw mechanism is operative: 

Tree level: m\ = m 2 = , 7713 ~ m 2 D /MR . (HI) 

That two masses are zero at the tree level is a consequence of Mo being a 1 x 3 matrix, 
which, therefore, maps two vectors onto 0; this feature is then operative in A4 d +m ()107|) 
and M v (fTUHJ) . 

At 1-loop level, m 2 becomes non-zero by neutral-scalar exchange, bit m\ remains zero: 

1 M 2 

1-loop: m! = , m 2 ~ — — m 3 -f In (M R /M ) 2 . (112) 

l07T z V z 

In this order-of- magnitude relation for mi |112j . the mass Mq ~ v is a typical scalar 
mass (there are three physical neutral scalars in this model). A general discussion for 
arbitrary numbers of left-handed lepton doublets, right-handed neutrinos singlets and 
Higgs doublets is found in |lllj . The full calculation of the dominant 1-loop corrections 
to the seesaw mechanism is presented in [113J. 

The model has the following properties. It predicts a hierarchical spectrum, therefore 
J Am|/ Am 2 tm ~ m 2 /m 3 e ~' 0.17. It gives the correct order of magnitude of 777-2/7713 
with 1/16tt 2 ~ 0.0063, M 2 /v 2 ~ 1, \n(M R /M ) ~ 10. The most interesting property of 
these 1-loop corrections to the seesaw mechanism is the fact that the suppression relative 
to the tree level terms is given solely by the loop integral factor 1/167T 2 |111[ I113j . The 
mixing angles are undetermined, but without fine-tuning they will be large in general. 
Consequently, the model has no argument for atm — 45° and small angle #13 ; these must 
be reproduced by tuning the parameters of the model. This is easily achieved because in 
good approximation | ZTq-3 | = \Mo a \l rn D |H2j . 

It is interesting to note that R-parity-violating supersymmetric models have a built-in 
seesaw mechanism where the heavy Majorana neutrinos are replaced by the neutralinos 
which are also Majorana particles. See, e.g., Ref. |114j and citations therein. 

4.5 A model for maximal atmospheric neutrino mixing 

Her we want to discuss a model based on tree-level seesaw masses. While in the previous 
model the emphasis was on explaining the ratio Arrig/Am^ and we assumed to obtain 
the mixing angles by tuning of model parameters, here we take the opposite attitude. 
One of the problems of obtaining maximal atmospheric neutrino mixing and large but 
non-maximal solar neutrino mixing by a symmetry is that the mixing matrix (jllOJl has a 
contribution also from the charged lepton sector. Now we introduce a framework where 
this problem is avoided by having Mr as the only source of neutrino mixing. 
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The framework: We start with the Lagrangian (|1U5|) and allow for an arbitrary number 
nff of Higgs doublets. Then how can one avoid flavour- changing interactions via tree-level 
scalar interactions? We assume that the family lepton numbers L e ^ tT are conserved in all 
terms of dimension 4 in the Lagrangian but L a is softly broken by the vr mass term. This 
allows one to kill several birds with one stone: Flavour-changing neutral interactions are 
forbidden at the tree level, Mi and Mr, are diagonal and neutrino mixing stems exclusively 
from Mr |115j . as announced above. However, soft breaking of lepton numbers occurs at 
the high scale nriR. It has been shown in Ref. |116j that this assumption yields a perfectly 
viable theory with interesting properties: All 1-loop flavour-changing vertices are finite 
because of soft L a breaking; for nu > 1, such vertices where the boson leg is a neutral 
scalar and the exchanged boson is a charged scalar do not decouple in the limit uir — > oo 
(for nil = 1 there is decoupling); diagrams with a 7 or a Z and box diagrams always 
decouple. As a consequence, the amplitudes of \i — > ej, Z — > e~ . . . are suppressed 
by l/m R , whereas, e.g., the amplitude of fi — > 3e tends to a constant for wir — > 00 and 
is suppressed — though much less than the previous amplitudes — because it contains a 
product of four Yukawa couplings. The decay rate of the latter process in this framework 
might be within reach of forthcoming experiments. For details see Ref. |llb'| . 



Maximal atmospheric neutrino mixing: Within the framework of soft L a -breaking 
we introduce now a symmetry: 

^2 : D^ L <-> D tL , v^r <-> u tR , fi R <-> t r . (113) 

This symmetry makes Mo and Mr ^-invariant and, therefore, transfers to the neutrino 
mass matrix M. v . Thus we obtain the result 




M v — \ y z w \ =^ U 



( cos9 sine \ 
sin e cos e 1 

V2 v% v% 

sing cos 9 1 

\ V2 V2 V2 ) 



(114) 



where 9 = On is the solar mixing angle. Summarizing, we have constructed a model where 
the solar mixing angle is free and without fine-tuning it will be large but non-maximal; 
atmospheric mixing is maximal; \U e ^\ = S13 = 0. 

The above results are stable under radiative corrections because the mass matrix ()114j) 
was realized by the symmetries U(1)l 6 x U(1)l x U(1)l t , which are softly broken, and 
I12, which has to be broken spontaneously in order to achieve m M 7^ m T . This can be 
done with a minimum of three Higgs doublets and an auxiliary Z' 2 , without destroying 
the form of the mass matrix (|114|) |115j . Since the Z 2 of Eq. ()113j) does not commute 
with L^t, the full group generated by our symmetries is non-abelian |117j . The essential 
features of this model can be embedded in an SU(5) GUT |117j . 



5 Conclusions 

In recent years we have witnessed great progress in neutrino physics. Eventually, it has 
been confirmed that the solar neutrino puzzle is solved by neutrino oscillations, first 
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conceived by Bruno Pontecorvo in 1957. At the same time, matter effects in neutrino 
oscillations — as occurring in the LMA MSW solution — have turned out to play a decisive 
role. It is firmly believed that also the atmospheric neutrino deficit problem is solved 
by neutrino oscillations though for the time being the final proof is still missing but is 
expected to be provided soon by LBL experiments. 

Despite the great progress, which provides us with a first glimpse beyond the SM, 
there are still many things we would like to know. For instance, Are neutrinos Dirac 
or Majorana particles? What are the absolute neutrino masses? Of what type is the 
neutrino mass spectrum? Do neutrinos have sizeable magnetic moments? 

As for field-theoretical models of neutrino masses and mixing, theorists are groping in 
the dark. As it has turned out, neutrino mass spectra and the neutrino mixing matrix 
are very different from the quark sector. Though there are many ideas, very few of them 
account naturally for some of the neutrino properties. Despite the big increase in our 
knowledge about neutrinos even basic questions for model building have no answer at the 
moment. Some of the basic questions are the following: Is the neutrino mass and mixing 
problem independent of the general fermion mass problem? Are neutrino masses small by 
radiative, seesaw or other mechanisms? Is the solution for the neutrino mass and mixing 
problem situated at the TeV scale or the GUT scale? Do we need a flavour symmetry for 
its solution and, if yes, of what type is this symmetry? 

We want to stress that it is no problem to accommodate neutrino masses and mixing 
in theories, but the problem is to explain the specific features for neutrinos. As for 
the neutrino nature, there is a theoretical bias toward Majorana nature, e.g. from the 
seesaw mechanism and GUTs, in particular, from GUTs based on 5*0(10). For the time 
being, there are simply not enough clues for a definite mechanism for neutrino masses 
and mixing. Among others, possible future clues provided by experiment would be an 
atmospheric mixing angle very close to 45°, which would point toward a non-abelian 
flavour symmetry; knowledge of the value of [U^l = S13 and the neutrino mass spectrum, 
which would teach us more about the mass matrix; discovery of scalars with masses < 
1 TeV, which would show that fermion masses are most probably generated by the Higgs 
mechanism; discovery of SUSY partners of ordinary particles, which would assure us that 
we have to take SUSY into account in model building. Evidence for flavour- changing 
decays like \i — > e^y, ji — > 3e, etc. would also give a valuable input. 

At any rate, ongoing and future experiments will continue to provide us with exciting 
results, which will enhance the prospects of constructing viable mechanisms for explaining 
the specific neutrino features. 
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